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Photon number resolving detectors can be highly useful for studying the statistics of multi-photon quantum 
states of light. In this work we study the counts statistics of different states of light measured on multiplexed 
on-off detectors. We put special emphasis on artificial nonclassical features of the statistics obtained. We show 
new ways to derive analytical formulas for counts statistics and their moments. Using our approach we are the 
first to derive statistics moments for multi-mode thermal states measured on multiplexed on-off detectors. We 
use them to determine empirical Mandel parameters and recently proposed subbinomial parameters suitable for 
tests of nonclassicality of the measured states. Additionally, we investigate subpoissonian and superbunching 
properties of the two-mode squeezed state measured on a pair of multiplexed detectors and we present results 
of the Fano factor and second-order correlation function for these states. 
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1. Introduction 

Multi-photon states of light are highly applicable in pre¬ 
cise quantum metrology (iH3. Among these states, the 
NOON states ii and the squeezed states [H are of the 
particular importance. The latter type of states have re¬ 
cently been successfully applied, for instance in increas¬ 
ing the precision of gravitational interferometer Q or in 
low-noise quantum imaging jH- 

Recent advances in technology allow for measure¬ 
ment of photon statistics of light with photon num¬ 
ber resolving (PNR) detectors. They have lived to see 
many implementations among which the most popu¬ 
lar are multiplexed on-off detectors based on the pho¬ 
ton chopping concept [H. They have been manufac¬ 
tured as fiber loop detectors |l0l4l2|. multi-pixel pho¬ 
ton counters (MPPC) (3 or as single photon-sensitive 
cameras BUI- Other types of detectors also in use 
include calorimetric transition-edge sensors 0 and hy¬ 
brid photo-detectors [iM3- The mutltiplexed on-off 
detectors have clear advantages such as a fast response 
and a relatively easy construction, since many of them 
are based on fast avalanche photodiodes. 

Up till now there have been a number of successful ex¬ 
periments using PNR detectors, thus developing further 
the knowledge on quantum states of light and its sources 
Many other experiments may be enhanced by 
the use of PNR detectors, for example, the observation 
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of macro-micro entanglement [22| or quantum imaging 

i- 

Here in this paper we will focus on the counting prop¬ 
erties of the multiplexed on-off detectors, which are the 
type of PNR detectors most often used. These detectors 
alter the counts statistics as compared to the photon 
statistics of light illuminating a detector [1, [23l - l25l |. A 
modification in counts statistics often leads to seemingly 
nonclassical properties of measured light [m,[HI; there¬ 
fore, an appropriate interpretation of the counts statis¬ 
tics is indispensable. 

In particular when it comes to an experiment, ad¬ 
equate criteria of nonclassicality based on empirical 
counts statistics have to be applied Criteria 

of nonclassicality for counts typically allow us to deter¬ 
mine qualitatively whether one is dealing with a quan¬ 
tum state or a classical one. In order to better differen¬ 
tiate between these states, the counts statistics models 
are indispensable. In many cases the measured state can 
be identified only based on the analysis of the mean and 
the variance of counts. 

In this paper we put a special emphasis on the mul¬ 
timode thermal states of light for which we are the first 
to derive the analytical model of counts statistics. Cur¬ 
rently there is a great interest in the community to ob¬ 
serve multimode light, not only generated in the sponta¬ 
neous parametric down-conversion, but also in atomic 
systems including room-temperature alkali metals va¬ 
pors (sol-l^. These experiments could be naturally 
reimplemented into PNR regime with the use of single- 
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photon sensitive cameras as in (33 |. 

Furthermore, in this article we present the numeri¬ 
cal model for two-mode squeezed states observed on two 
identical PNR detectors. We present the results for non- 
classical measures such as the Fano factor and the sec¬ 
ond order correlation function, commonly used in exper¬ 
iments. 

2. Detector counting properties 

Multiplexed detectors consist of a finite number of N 
Geiger-type on-off detectors. Number of detectors N 
is the parameter which determines the counting statis¬ 
tics of the detector. N is also the maximum number of 
counts that can be measured. 

The detector operation can be described by giving the 
conditional probability of obtaining k counts provided 
that the detector was illuminated by n photons 

PN{k\n) = ^ k\S{n,k), (1) 

where S{n, k) = -g X]i=o(—1)* 

ling number of the second kind j^. For instance, Eq. 
m yields PAr(0|n) = 0 for all n > 1 and p{l\n) = • 

Typically it is enough to know the conditional proba¬ 
bilities pN{k\n) for a finite value of n since the detector 
saturates for large n as we illustrated in Fig. [T] (a-b). 

Due to finite quantum efficiency r] of the detector the 
model Eq. [T] has to include photon losses p^(n|m) = 

^n(i 

This can also be done analytically for specific n and 
k, for instance Pn{^\'^) = (1 “ 

»(ii")=^^(1 V )"-!)- (2) 

We have gathered a whole set of conditional probabil¬ 
ities in Eig. [1] (c-d). 

Here in this paper we derive an analytical and numeri¬ 
cal model for an idealized detector subject to certain as¬ 
sumptions. We shall assume that the component detec¬ 
tors are very similar, thus principally they will have the 
same quantum efficiency 77 . Moreover, the dark counts 
and the cross talk between the component on-off detec¬ 
tors will be excluded from our model. 

If one cannot apply these assumptions to the used 
detector, another, alternative approach should be used 
such as Bayessian approach to treat the cross-talk [TsI, 
[ 3 ^, Including all imperfections in a general form to 
the theoretical model would make it very complicated 
and of little use. Instead, one can find the values of 
PN{k\n) through detector tomography (^. [^.[38l| . 

3. Counts statistics on a single detector 

In the experiment with the PNR detector, we have 
access to the history of counts collected by the de¬ 
tector, which yields empirical probability distributions 


/ m 




Figure 1. Values of conditional probabilities pN{k\n) of get¬ 
ting k counts when a detector is illuminated by exactly n 
photons, for a detector consisting of N multiplexed compo¬ 
nent on-off detectors with quantum efficiency 77 . (a) V = 2 , 
77 = 1 (b) (a) iV = 4, 77 = 1 (c) V = 4, 77 = 0.8 (d) V = 4, 

T] = 0 . 2 . 


for counts statistics c/^. It can be expressed as Ck = 
where fn is the photon number distri¬ 
bution in the measured state. 

Detector losses can be taken into account both by an 
appropriate modification of the conditional probabilities 
of the detector pN{k\n) and, as described earlier, in the 
photons statistics. 

In the following discussion we shall assume the statis¬ 
tics fn after losses and the detector will be described by 
no-loss conditional probabilities as in Eq. O- 

In the experiment we often confine ourselves to track¬ 
ing only the mean and the variance of the counts statis¬ 
tics. Determination of these values is often sufficient to 
identify the measured state. Now, we shall derive the 
moments of probability distributions for the counts. 

In general we can find average values of polynomial 
functions x(/c) of counts: 

{x{k)) = '^x{k)iAk\^—S{n,k)fn (3) 

/c=0 ^ ^ n=k 

To evaluate Eq. we can use two useful properties 
of Stirling number of the second kind (35l |: 

1 


J2S{n,k)x^ = -^^ (5) 

n=k 
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<k> Qp 


- <(Ak)^> - Qm=Qb=0 



(a) (b) 


which can be readily associated with the single mode 
thermal state of the mean (n) for which a = (1 + 
and b = (n)/(l + (n)). 

Then, using the property Eq. m, we find the average 
of any function of counts for given statistics fn- 

ix{k))f=a^x{k) 


Figure 2. (a) Mean and variance of counts for coherent state We focus particularly on the first and second moment 

on a detector with N = 4. (b) Measured Qf parameter is of the counts statistics; 
artificially nonclassical contrary to Qb = Qm = 0. 


where (•)n denotes the Pochhammer symbol: {y)n = 
T{y + n)/T{y). 

These two properties facilitate the determination of 
means and variances for states of light frequently used 
in experiments. For instance, using the property Eq. 61 
one can find the mean {k) and the variance {{Ak)‘^) of 
counts for the coherent state 

j n 77,1 


(k) = N{1 - (6) 

{(Akf) = N{1 - (7) 


{k)f 


obN 

{h-l){b{N -1) - N) 


( 10 ) 


{k^)f = 


ab{b + 1)N‘^ 


{b - l){b{N - 2) - N){b + N-bN)' 


( 11 ) 


We are also able to find moments of counts distribu¬ 
tion for another type of statistics yielded from Eq. ([8]), 
in particular: 


9n 


n\ 


which can be expressed as: 


( 12 ) 


Here we see that the detector reduces both the mean 
and the variance of counts as compared with the values 
for the photon statistics Fig. [2] (a). 

This reduction leads to the seemingly nonclassical 
properties of the measured light. One of nonclassical- 
ity criteria for single-mode light is the negativity of 
the Mandel parameter Qm = ((An)^)/(n) — 1. The 
parameter can also be evaluated for counts statistics 
Qf = {{Ak)‘^)/{k) — 1, here for coherent states being 
always negative Qf = — 1 < 0 (Fig. [21(b)). 

A recently proposed modified criterion of nonclassi- 
cality is based on the subbinomial parameter mill: 


Qn = adb^mb'^^'^ = db,mb”^fn- (13) 

Note that Qn can be readily associated with the multi- 
mode thermal state with an average number of photons 

(n) and the number of modes M. is described by the 

\ M 

statistics gn Eq. (US! for a = ( i^J^m ) ’ ^ = 

{n)/{A4 + (n)) and m = M. — 1 (3^ . 

To calculate the moments for given statistics we 
simply do the following: 


Qb=N 


{k){N-{k)) 


- 1 , 


negative only for nonclassical states and for coherent 
states reconstructing the true value of Qm Mandel pa¬ 
rameter (Fig. [2](b)). 

4. Multi-mode thermal states of light 

Now we proceed to calculations for single and multi- 
mode thermal states of light which can viewed as a single 
subsystem of the squeezed states of light. 

At first, let us consider the following photons statis¬ 
tics: 


(8) 


{x{k))g = db,mb^{x{k))f (14) 

In this way, we can find the mean and the variance for 
a single-mode thermal state, wherea = (1 + and 

b = (n)/(l + (n)), similarly as in (^ : 


(^)Th,l 


{n)N 
(n) + A 


((A/c)^)Th,l 


{n)N‘^{{n)N^{n)^N) 
((n)+7V)2(2(n) + 7V) 


Equations da, (nni) and m and substitutions for a 
and b yield analytical expressions for the mean {k)Th,M 
and the variance {{Ak)‘^)Th,M foi* ^^ch M. and N. For 


fn = ab^, 


















4 



(a) (b) 


Figure 3. Theoretical values for empirical Qf parameters 
for (a) a single-mode thermal state, (b) multi-mode thermal 
state A4 = 4 plotted versus mean number of photons (n) for 
different N. Depending on the parameters of the state and 
the detector counts statistics may seem classical or artificially 
nonclassical. 


- <k> - <(Ak)2> 


- Qf - Qb - Qm 



<n> <n> 

(a) (b) 


Figure 4. (a) Mean and variance of counts, (b) input Qm, 
measured Qf and binomial Qb parameters versus mean 
number of photons in two-mode {M — 2) thermal state on a 
detector with N = 8. 


example, we give analytical expressions for {k)Th,M foi* 
two-mode A4 = 2 thermal state: 


(^)Th,At=2 = 


(n)N((n)-h4N) 

((n)+2Ar)2 


Further the second moment (k^}Th,M =2 can be ex¬ 
pressed analytically: 


(^^)Th,At =2 = 

_ (n}N^((n)^+6(n}^N+3(n)N(2N+l)+4N^) 

~ ((n)+Ny^((n}+2NP 

Analytical formulas for a higher number of modes can 
be readily found using any symbolic computation soft¬ 
ware and instead of presenting directly the formulas we 
gather the results for higher number of modes on plots. 

Having found means and variances parameters, we 
can construct Mandel parameters for counts and com¬ 
pare them with theoretical values Qm = (n}/A4. In 
Fig. [3] we show the empirical Qf parameter for detec¬ 
tors of different N for single- and multi-mode thermal 
state (A4 = 4). Depending on the mean number of in¬ 
put photons (n) and N, the measured Qf parameters 
appear to be nonclassical is certain regimes. 


5. Counts properties of squeezed states 

The above results can be measured in a single subsystem 
for single- and multi-mode squeezed states using a single 
detector. If we have two detectors, we can calculate the 
joint counts statistics Ck^^k 2 related to the input statistics 
/ni,n 2 by relation: 


oo oo 

Cki,k2 = E E Fiv(fci|ni)pjv(fc2|n2)/ni ,n2 

ni=ki n 2 =k 2 

Calculations of the joint count statistics allow us to 
detect subpoissonian correlations between the two sub¬ 
systems. A good measure to quantify these correlations 
is the Fano factor [T3| which, evaluated for counts, can 
be expressed as: 


_ {{Ajk.-hm 

ih) + {k2) ’ 

which is R> 1 for all classical states. 

Another commonly used parameter to characterize in 
experiments the two modes states of light is the second 
order correlation function, which evaluated on counts 
statistics can be expressed as: 

( 2 ) _ (^ 1 ^ 2 ) 

^ {ki){k2)’ 

It can be used to test if the investigated light demon¬ 
strate the superbunching properties i.e. for > 2 , 
which is characteristic for squeezed vacuum states (40| 
but in principle may be also obtained using interference 
of thermal states . 

Now we shall focus on properties of joint counts statis¬ 
tics of the two-mode squeezed state of light: 


00 

I'*/’) = Vi - icp 

n=0 

where is the squeezing parameter. The Fano factor 
for such a state without losses equals i? = 0 , whereas 
the second orderd correlation function becomes always 
over two: = 1 + 1 /|CP ^ 2 . 

In Fig. Owe show the effect of a finite number N of 
component detectors and influence of losses on the joint 
counts statistics Ck^^k 2 i C = 0 - 8 - 

Both the finite number of detectors and the losses have 
influence on reducing the correlation between counts 
measured on two detectors. To better understand this 
mechanism, we found the Fano factors for different N^ 77 
and Q Fig. [ 6 ] (a) presents pure influence of a finite N^ 
without the losses. 

For small squeezing parameters C high number of 
component detectors N ensures a smaller Fano factor, 
whereas for high the effect is opposite. This is the ef¬ 
fect of artificial increasing of correlations due to the de¬ 
tector saturation. The saturation effects are even more 
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n, n, 


Figure 5. Two mode squeezed state with ( = 0.8 measured by 
two multiplexed detectors, (a) = oo, 77 = 1 . (b) = 00 , 

Tj = 0.8. (a) A/" = 4, ?7 = 1. (a) A/" = 4, 77 = 0.5. 




Figure 6 . Fano factors versus squeezing parameter for dif¬ 
ferent numbers of N detectors with (a) 77 = 1 (b) 77 = 0.5. 
A lower number of detectors N can artificially increase the 
correlation and consequently decrease the Fano factor. 


significant for the detector with finite quantum efficiency 
(Fig. [n](b)), where 77 = 0.5. Here only a smaller number 
of detectors N has an effect on reducing the Fano factor 
for each^. 

In all these cases, the calculated counts statistics al¬ 
ways remain nonclassical. In order to perform a reliable 
test of nonclassicality with no a priori knowledge of the 
state of light, one can apply the recently proposed crite¬ 
ria ( 2 §|. 


It is also instructive to view how the second corre¬ 
lation function is modified due to the limitations 
introduced by a finite number of component on-off de¬ 
tectors N. It can be seen in Fig. [7] that the low num¬ 
ber of component on-off detectors N decreases by 
no more than 1. This means that limit between super- 
bunched states and the bunched states = 2) will be 
exceeded for a sufficiently low N for squeezed states of 
the squeezing parameter higher than \(\‘^ > 1/2 as it is 
exemplified in Fig. [7](b). 



N N 

(a) (b) 


Figure 7. Second order correlation function g^‘^^ versus dif¬ 
ferent numbers N of detectors for 77 = 1 for \Q = 0.2 (a) and 
for |C| = 0.8 (b) for two-mode squeezed states. Finite num¬ 
ber N of component on-off detectors decreases g^‘^^ which for 
two-mode squeezed state always is higher than two, g^‘^^ > 2 . 
g^^^ is decreased by no more than one. For > 1/2 and 
sufficiently low N, the observed counts statistics does not 
preserve its initial superbunching properties as in (b). 


6. Summary 

In this paper, we investigate the statistics for counts 
measured by a detector illuminated by squeezed or ther¬ 
mal states of light, the latter being a traced subsystem 
of the squeezed state. 

In particular, we provide a universal manner of deriv¬ 
ing statistics moments which we have used to determine 
analytical formulas for multi-mode thermal states. In 
this work we put special emphasis on the issues related 
to artificial nonclassicality of measured statistics in the 
context of Mandel and subbinomial parameters. 

On the other hand, we show the influence of increas¬ 
ing and decreasing subpoissonian counts correlations for 
two-mode squeezed states of light measured on two mul¬ 
tiplexed detectors. We show that the Fano factor al¬ 
ways remains nonclassical. Moreover we provide the re¬ 
sults for the second-order correlation function g^‘^^ which 
could drop below two in certain cases. 

The results of the work can be used to identify the 
states of light measured by multiplexed on-off detectors, 
in particular to determine the number of modes in ther¬ 
mal states. The results can significantly contribute to 
designing experiments using multiplexed on-off detec¬ 
tors. They will play an increasingly important role in 
the future since the state of the art, photon number re¬ 
solving detectors are becoming a very efficient tool for 
studying the multimode, multi-photon quantum states 
of light. 
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